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Abstract 

Nuclear pairing gaps of normally deformed and superdeformed nuclei are investigated using the 
particle-number conserving (PNC) formalism for the cranked shell model, in which the blocking 
effects are treated exactly. Both rotational frequency w-dependence and seniority (number of 
unpaired particles) z^-dependence of the pairing gap A are investigated. For the ground-state 
bands of even-even nuclei, PNC calculations show that in general A decreases with increasing uj, 
but the w-dependence is much weaker than that calculated by the number-projected Hartree-Fock- 
Bogolyubov approach. For the multiquasiparticle bands (seniority u > 2), the pairing gaps keep 
almost w-independent. As a function of the seniority v, the bandhead pairing gaps A(z^, a; = 0) 
decrease slowly with increasing u. Even for the highest seniority u bands identified so far, iS.{v,oj = 
0) remains greater than 70% of A(i/ = 0,a; = 0). 
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I. INTRODUCTION 



Since the seminal article by Bohr, Mottelson, and Pines significant effects of nu- 
clear pairing were established in fundamental nuclear properties 0. Soon afterwards, the 
Bardeen-Cooper-Schrieffer (BCS) theory for metallic superconductivity and quasiparticle 
(qp) formalism were transplanted in nuclear structure literature to treat nuclear pairing 
correlation Now the BCS or more elaborate Hartree-Fock-Bogolyubov (HFB) ap- 

proximations are the standard methods in nuclear physics. However, along with their great 
successes, both BCS and HFB approximations for nuclear pairing raise some concerns ||^, |^. 
One of them is the non-conservation of the particle-number. Because the number of nucleons 
in a nucleus is not very large (n ~ 10^), particularly the number of valence nucleons (n ~ 10) 
dominating the nuclear low-lying excited states is very limited, the relative particle-number 
fiuctuation, 5n/n, is not negligible. Indeed, it was found that in all self-consistent solutions 
to the cranked HFB equation a pairing collapsing occurs for angular momentum I greater 
than a critical value Ic [§. 

To restore this broken symmetry, many works have been done. The Lipkin-Nogami (LN) 
method [pHTlH has been quite extensively used in these years. After using this approximate 
particle-number projection method, pairing phase transition disappears [|12], |13[. However, 



earlier studies showed that the LN method broke down in the weak pairing limit [0, |I5|. At 
the same time, various particle-number projection approaches of pairing interaction in BCS 



or HFB formalism have been developed |p!6|-p!9|. In these approaches, the ideal treatment is 
variation after projection, but when spin goes higher, this method becomes very complicated 
and computational expensive. All these methods tried to solve the problem of the particle- 
number nonconservation. However, when it is achieved, no pairing phase transition was 
found [jl6|, [l^. This proves that the occurrence of nuclear pairing collapsing originates 
from particle-number non-conservation. Other than the variational approach, some methods 
directly solving the corresponding Schodinger equation have been developed P, |20| . In these 
methods, the particle-number is strictly conserved. The particle-number conserving method 
used in Ref. pO|, in which the single-particle states stem from the Hartree-Fock mean field, 
is a little different from the method used in our work, in which the single-particle states 
stem from the Nilsson model. 

Another problem related to the violation of particle-number conservation is the occurrence 
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of spurious states in the BCS (HFB) qp formalism. As pointed out by Richardson an 
important class of low-lying excitations in nuclei cannot be described in the standard BCS- 
or HFB-like theories. The remedy in terms of the particle-number projection considerably 
complicates the algorithm, yet failed to properly describe the higher energy spectrum of the 
pairing Hamiltonian |^. 

The most concerned issue is the proper treatment of the Pauli blocking effect on pair- 
ing, which is responsible for the odd-even differences in nuclear properties (binding energies, 
moments of inertia, etc.) As emphasized by Rowe [^], while the blocking effects are straight- 
forward, it is very difficult to treat consistently in the qp formalism because they introduce 
different qp bases for different blocked orbitals. Indeed, it was shown that the properties of 
a rotational band are very sensitive to the Coriolis response of the blocked single-particle 
orbitals [2^ . 



In this article, to investigate the pairing reduction due to rotation and blocking we use the 
particle-number conserving (PNC) formalism for treating the cranked shell model (CSM) 
with pairing interaction, in which the particle-number is conserved and the blocking effects 
are treated exactly. The details of the PNC formalism for calculating the MOI have been 
given in [Q. Only the PNC formalism for calculating the nuclear pairing gap is given 
in Sect. II. Sect. Ill gives PNC calculations for nuclear pairing gaps of various types of 
pair-broken rotational bands in normally deformed (ND) and superdeformed (SD) nuclei 
(seniority (number of unpaired particles) z/ > 2 for even-even nuclei, u > 1 for odd-A 
nuclei), as well as the rotational frequency cj- dependence and seniority i^-dependence of 
pairing gaps. A brief summary is given in Sect. IV. 

II. PNC FORMALISM FOR NUCLEAR PAIRING GAP 

The CSM Hamiltonian of an axially deformed nucleus in the rotating frame is 

-f^csM = Ho + Hp = Hj^ii — uJx + Hp , (1) 
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where i^Nii is the Nilsson Hamiltonian, —uJx is the Coriohs interaction with cranking fre- 
quency u! about the x axis, Hp = Hp{0) + Hp{2) is the pairing interaction, 

Hp{0) = -Go^a^a^afja^ = -Gq^s^s^ , (2) 
Hp{2) = -G2^q2{0(l2{v)at"'t^f^^v = -G2^q2{0(l2{v)spn ^ (3) 

where ^ (f/) labels the time-reversed state of a Nilsson state ^ (77), q2iO = A/l67r/5(^|r^y2o|0 
is the diagonal element of the stretched quadrupole operator, and Gq and G2 are the effective 
strengths of monopole and quadrupole pairing interactions respectively, = a^a^ {sr^ = 
afjarj) is the pair creation (annihilation) operator. 

In the PNC calculation, Hqsm is diagonalized in a sufficiently large cranked many-particle 
configuration (CMPC) space |^ and is expressed as 

l^) = $^aK) (a real), (4) 

i 

where |i) is an eigenstate of Ho with configuration energy eI^\ characterized by the particle- 
number N, parity vr, signature r (= e~*'^") and seniority u (number of unpaired particles). 
For the seniority u = ground state band {K"^ = 0''") of an even-even nucleus (qp vacuum 
in the BCS formalism), each \i) in Eq. |^ is of the product form 

44---I0), ^^r^^.... (5) 

For the seniority u = 1 band (~l-qp band in the BCS formalism) in an odd-even nucleus, 
\i) is of the form 

«^44---|0), e7^^^---(^A) , (6) 

where A is the blocked single-particle state, i ^ f] ^ ■ ■ • \) (Pauli blocking effect) and the 
angular momentum projection along nuclear symmetry z-axis K = Qx. For the seniority 
u = 2 band (~2-qp band in the BCS formalism) in an even-even nucleus, \i) is of the form 

4«^4<---|0) ' , (7) 

where X ^ a are two blocked single-particle states {K = Qx + ^a)- The PNC forms of 
the u > 2 (multiquasiparticle) bands are similar. Strictly speaking, due to the Coriolis 



interaction —uJx, v and K are not exactly conserved for a; 7^ 0. Walker and Dracoulis ||25 



pointed out that some forms of ii'-mixing must exist to enable the /T-forbidden transition 
observed in a lot of low-lying rotational bands of axially symmetric nuclei. However, in the 
low-w region, v and K may be served as useful quantum numbers characterizing a low-lying 
excited rotational band. 



The kinematic and dynamic MOIs for the state |\E') follows pi, |2| 



J« = -(vl>| J.|vl>) , J(2) = A(vl/| J,|vp) , (8) 

where 

(*|J,|^) = ^Cf(z|J,|z) + 2^aC,-(z|J,|j) , (9) 

i i<j 

is the angular momentum alignment of the state |\E'). 



In the PNC formalism, the nuclear pairing gap may be reasonably defined as [|16|, 0, |27 

1 



Go 



(10) 



where l^') is a PNC eigenstate (Eq. ^ of Hqsm with eigenvalue E. In the BCS formalism 
for HcsM with the monopole pairing interaction only. Hp = —GoS'^S, where S'^ = J2^^t^ 
S = J2r] ^^'^ 'i'P vacuum band |0)) 

\0)) = U^{U^ + V^s+)\0) , Ul + Vl = l, (11) 
A is reduced to the usual definition of nuclear pairing gap A 

A = GomS+m=GoJ2u^V^- (12) 

Calculations show that for the low-lying excited eigenstates of Hqsm, the number of 
important CMPC's (with weight > 1%, say) is very limited (usually < 20 for the ND rare- 
earth nuclei), thus it is not difficult to get sufficiently accurate solutions to the low-lying 
excited eigenstates of Hqsm by diagonizing i^csM in a sufficiently large CMPC space 
To ensure the PNC calculations for nuclear low-lying excited states both workable and 
accurate 0, |2^, it is essential to adopt a CMPC truncation (Fock-space truncation) in the 
PNC calculation in place of the usual single-particle level (SPL) truncation in shell model 
calculations. This is understandable from the perturbation expansion of Hqsm (H); as it refers 
to a many-particle system with pairing interactions. In general, the lower the configuration 
energy of the MPC is, the larger the weight of the corresponding MPC in low-lying excited 



eigenstates of Hqsm will be. The stability of the final results with respect to the basis cut-off 



has been illustrated in details by Molique and Dudek as well as in p9[| . 

In the following PNC calculations, Hqsm is diagonalized in the CMPC space with dimen- 
sion 1500 for both protons and neutrons. The corresponding effective proton and neutron 
pairing strength are adopted to reproduce the experimental odd-even differences in nuclear 
binding energies. Proper Nilsson level schemes are adopted to reproduce the experimental 
bandhead energies and MOI of the low-lying excited seniority u = 1 (1-qp) bands. Thus, 
the pairing gaps A of various low-lying excited bands can be convincingly extracted by the 
PNC calculations without any free parameter. 

III. CALCULATIONS AND DISCUSSIONS 

In this section the PNC calculations for nuclear pairing gaps of some typical bands in 
ND and SD nuclei are presented. The rotational frequency w-dependence and seniority 
z/-dependence of the pairing gaps are discussed in detail. 

A. Ground state bands of ^^^Yh and ^^SHf 

The angular momentum dependence of pairing gaps A„ (neutrons) and Ap (protons) for 
the ground state band (gsb) of i^syb i^sjjf have been calculated in the number-projected 
HFB (NHFB) approach in Ref. Jl^. The pairing gap reductions in the observed angular 
momentum range 1 = 0^ AAk (^ w = 0.61 MeV/h) for ^'^'^Yb(gsb) and 1 = 0-)- Mh 
{^uj = 0.52 MeV/n) for i68Hf(gsb) calculated by NHFB are [|16|, |l7| 

^^(^ = ^^^^^63%, for-Yb(gsb) 
A„(/ = 0) Ap(/ = 0) 

K{I = Uh) ^^g^^^ for ^^«Hf(gsb) (13) 



A„(/ = 0) 

For comparison, in this section the w-dependence of pairing gaps of ^^^Yb(gsb) and 
^^^Hf(gsb) are calculated using the PNC formalism. To validate the PNC calculations of 
Ap (proton) and A„, (neutron) (Fig. 0(b)), the kinematic MOIs J^^-* are also calculated un- 
der the PNC formalism and compared with the experiments I^D], ^ (see Fig. |T](a)). The 



experimental MOIs J'-^^ are very well reproduced by the PNC calculations (except in the 
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FIG. 1: The MOIs and pairing gaps A for the ground state bands of ^^^Yb and ^^^Hf. (a) The 



experimental MOIs [30, 31| are denoted by the soUd circle •, and the calculated MOIs by the PNC 
method are denoted by solid lines. The Nilsson parameters (k, /i) and deformation {e2,£i) are 
taken from p^ , |3^ . The monopole and quadrupole pairing strengths for protons and neutrons are 
adopted to reproduce the odd-even differences in nuclear binding energies, Gn = 0.30, G2n = 0.010, 
Gp = 0.29 for i68Yb; G„ = 0.39, Gp = 0.35 for ^^^Hf. (b) The PNC calculated pairing gaps for 
protons (neutrons) are denoted by solid (dashed) lines. 

bandcrossing region). Thus, we believe the PNC calculations of pairing gaps (w-dependence, 
z/-dependence, etc.) are trustworthy. In the observed rotational frequency range, the pairing 
gap reductions calculated in PNC formalism are 

O.eiMeV/^) _ Ap{uj = 0.61MeV/;i) 



A„(w 



A„(a; = 0) 
A„(a; = 0.52MeV/;i) 



70%, 
^ 70%, 



Ap{uj 



Ap{uj 



0) 

0.52MeV/;i) 



for ^^^Yb 
83%, for i^^Hf 



(14) 



0) ' Ap{uj = 0) 

which remains more than 70% of the bandhead value in all experimental u range. As 
expected, in both NHFB and PNC formalism no pairing phase transition from superfluidity 
to normal motion (A — 0) is found with increasing u. However, the w- dependence of 
A in PNC calculations is weaker than that calculated by NHFB approach. By the way, 
it is noted that due to the neutron sub-shell effect at = 98, in both PNC and NHFB 
calculations for ^^®Yb, the pairing gap reduction of neutron is larger than that of proton. It 



was noted by I. Hamamoto [34 that an inherent issue of CSM is the violation of rotational 



symmetry, and the reliability of calculations in the CSM, particularly in the bandcrossing 
region, is questionable. Afterwards, the angular momentum projection techniques were 
developed |^5[. It is interesting to note that the w-dependence of the pairing gaps for the 
gsb of ^^^Yb calculated by the angular momentum projection technique are similar to 
that of the PNC calculations. 



B. Multiquasiparticle bands of the heavier rare-earth nuclei {A ~ 178) 

The seniority i^-dependence of nuclear pairing gaps have been investigated by Dracoulis 
et al., using the LN method [0. They showed that the bandhead pairing gap A(z/, a; = 0) 



decreases approximately by 

A{u, 00 = 0) = {0.75y^^A{u = 0, w = 0). (15) 

In this section we will investigate the i^-dependence of the A using the PNC formalism. To 
get pairing gaps for these multiquasiparticle bands, a proper Nilsson level scheme for the 
deformed heavier rare-earth nuclei {A ~ 178) is necessary. However, the level scheme (Lund 



systematics) |33| is unable to properly reproduce the experimental bandhead energies of 
the low- lying excited 1-qp bands of ^^''Ta, particularly the gsb, 7r7/2"''[404]. So the Nilsson 
parameters (k, //) in are slightly adjusted (see the caption of Fig. ||). Fig. || shows the 
experimental ||38| and calculated MOIs of four 1-qp bands in ^^^Ta. 

In Fig. ^ we show the PNC calculations of the proton pairing gaps for i/ = 1,3 bands 
in ^^''Ta and the gsb, and multiquasiparticle bands with Up = 2 and 4 proton configurations 
in i^^W fS^-il. The cu-dependence of Ap for ^■^®W(gsb) is similar to that of i^*Yb(gsb) 
and ^^*Hf(gsb) (see Fig. [l|(b)). The cj-dependence of Ap's of the four low-lying excited 1- 
quasiproton bands of ^'^'^Ta is similar to ^^^W(gsb), so only the Ap for the gsb (7r7/2"'"[404]) 
of ^'^''Ta is shown in Fig. ^. The dependence of the z/p = 2 configuration ir^S^ (7/2"*" [404] (8) 
9/2^ [5 14]) in ^^^W is weaker. For the low- lying excited multiquasiparticle (seniority Up > 2) 
bands, Ap's keep nearly u -independent. In fact, for realistic nuclei, the blocking effects on 
pairing are significant only for a few orbitals nearest to the Fermi surface. For low-lying 
excited multiquasiparticle bands, a few orbitals nearest to the Fermi surface are almost 
blocked, and for orbitals far from the Fermi surface the w- dependence of the blocking effects 
on pairing are quite small. 
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FIG. 2: MOIs of four low-lying seniority v = \ bands in ^^''Ta. The experimental MOIs [38| are 
denoted by • (a = 1/2) and o (a = —1/2), respectively. The calculated MOIs by the PNC method 
are denoted by solid lines (a = 1/2) and dotted lines [a = —1/2), respectively. The Nilsson 
parameters (k, ^) in are slightly adjusted to reproduce the bandhead energies of the 1-qp 
bands. For protons, K4 = 0.060 {N = 4), K5 = 0.061 {N = 5), fi^ = 0.55, ^5 = 0.69. For neutrons, 
K5 = 0.066, kq = 0.058, /is = 0.49, fiQ = 0.40. The deformation parameters (£2, e4)=(0.24, 0.04) 



are from [33|, i.e., an average of the neighboring even-even Hf and W isotopes. The effective 
pairing interaction strengths for both proton and neutron, Gn = 0.26 MeV, Gp = 0.26 MeV, are 
determined by the experimental odd-even differences in nuclear binding energies. 

As a function of seniority z/, in general, the pairing gap A(z/) gradually decreases with 
increasing u. The pairing gap reductions at bandhead {u = 0) calculated by PNC method 
are 



Ap(z/ = 1) 



0) 
= 3) 
0) 



91% , 
i80% 



Ap(z/ = 2) 
Ap{iy = 0) 
Ap(z/ = 4) 
Ap{u = 0) 



86% , 
78% . 



(16) 



|37[] . Even for the highest seniority u bands 



which are weaker than that given in Eq. (^ 
identified so far, the pairing gaps A{iy,u = 0) is always larger than 70% of A(z/ = 0, w = 0). 
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FIG. 3: The protons pairing gaps Ap for the fp = 1, 3 bands in Ta (dotted hnes) and Up = 0, 2, 4 
configurations in ^^^W (solid hnes).^'^^Ta: Up = 1 band (gsb), 7r7/2+[404]; i/p = 3, K'^ = 17/2+ 
band at 1523 keV ||], 7r3l7/2+(7/2+[404] (g) 9/2-[514] ® l/2-[541]). ^^^^W: gsb, Up = 0, K"" = 0+; 
i/p = 2 configuration 7r28~(7/2+[404] (g) 9/2" [514]) in i^'^ = 15+ (vr^S" 1/^7-) band at 3653 
keV; i/p = 4 configuration 7r4ll-(7/2+[404] 5/2+[402] (g 9/2-[514] ® l/2-[541]) in K'^ = 18+ 
{ttHI- (g) i^^7~) band at 4878 keV |3^, |9|. 

C. Ground state bands of 238u and ^^^No 

The PNC calculations for the gsb's of actinide nuclei, ^ss-q ^-^^^ '^^^'Mo, are shown in Fig. ^. 



The experimental data of MOIs are taken from The PNC calculations for the pairing 

gap reduction show 



A„(a; = 0.30MeV/n) Ap{uj = 0.30MeV/n) 

- ~ 85% , - 



An{u = 0) 

A„(a; = 0.30MeV/n) 

An{uJ = 0) 



Ap{u = 0) 



for ^•'^U(gsb) 



m. ^^(^^H^MfXM « 76% , for -No(gsb) (17) 

Ap{uj = 0) 



i.e., Ap's and A„'s decrease very slowly with increasing u, quite similar to the rare-earth 
nuclei. 
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FIG. 4: The MOIs and pairing gaps of the gsb of ^^^U and ^^^No. (a) The experimental MOIs ||4^, 
41 1 are denoted by sohd circle • (a = 0, 1/2) and open circle o (a = —1/2). The calculated MOIs 
by the PNC method are denoted by solid lines (a = 0,1/2) and dotted lines {a = —1/2). The 
Nilsson parameters {k, fi) are taken from |32] both for ^38^ ^^^^ ^^^No. Deformations (£2, £4) = 
(0.21,-0.04) of are taken from |33|. For ^^^No, the quadrupole deformation parameters for 
the neighbor even-even nuclei are deduced as P2 = 0.28 it 0.02 from the experiment Here we 
choose £2 = 0.26 and and £4 = 0.01 in our calculation for ^^^No. To reproduce the experimental 
odd-even difference in binding energies and the w-dependence of MOIs, it seems that a small 
amount of quadrupole pairing interaction is necessary, i.e., G„ = 0.29 MeV, G2n = 0.015 MeV, 
Gp = 0.32 MeV, Gap = 0.040 MeV for 238u; G„ = 0.22 MeV, G2„, = 0.010 MeV, = 0.26 MeV, 



G2p = 0.010 MeV for ^^-^No. (b) The PNC calculated pairing gaps Ap's (A^'s) for the gsb's of 
^^®U and ^^'^No are denoted by solid (dashed) lines. 

D. Pairing gaps of SD bands 



Chasman |Q pointed out that BCS treatment of nuclear pairing is not appropriate for 
SD states because the single-particle level density near the Fermi surface is low and the BCS 
method is not correct in this limit. In most cases, the linkage between the SD bands and 
low-lying excited states have not yet been established, thus the actual spins of some SD 
bands are not determined, but the dynamic MOIs J^'^^ can be extracted from the observed 
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FIG. 5: The pairing gaps of SD bands in Hg isotopes calculated by the PNC formalism. The Nilsson 
level scheme are taken from [43]. "'^^^Hg(l): 1-quasineutron SD band, z^5/2^ [512], a = 



-1/2; 



i93Hg(2b): 1-quasineutron SD band, z^9/2+[624], a = 1/2. i*^^Hg(l): quasi-vacuum SD band 
(a = 0); i94Hg(2): 2-quasineutron SD band i/27-(5/2- [512] ® 9/2+ [624]), a = 0. 



differences in subsequent 7 transition energies, 

/2)(J) = Aff/ [E^{I + 2^1)- E^{I ^1-2)] . 



(18) 



Obviously, the accuracy of J(2) is lower than J^^\l) = {21 + 1)^7 [E^{I + 1 ^ / - 1)]. 
However, the actual spins of some SD bands have been established experimentally (e.g., see 
Ref. m), thus the J^^^'s can be accurately extracted. The cu-dependence of experimental 
MOIs for series of SD bands l|46|-48l were reproduced very well by the PNC calculations 

In 



50 



for the CSM with both monopole and quadrupole pairing interactions [^, | 
Fig. ^, the PNC calculations of the pairing gaps A for the z/„ = 1 SD bands ^^^Hg(l) 
(z/5/2-[512],a = -1/2), i93Hg(2b) (z/9/2+[624], a = 1/2), the u = SD band ^^^Hg(l), and 
= 2 SD band ^^^Hg(2) (z/^T", 5/2-[512] ® 9/2+[624],a = 0) are shown in Fig. |. The 
Nilsson level scheme are taken from [Q. PNC calculations show that: 

(a) For SD bands in Hg isotopes, Ap (proton)^ A„ (neutron), which is caused by the large 
gap at Z = 80 in the proton Nilsson level scheme of SD Hg isotopes. 

(b) For SD bands, no pairing collapsing is found with increasing co either. For ^^^Hg(l) 
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(z/ = 0, SD band), the pairing gap reduction with increasing u is, 

A„(a; = 0.50MeV/;i) ^^^^ A,(a; = 0.50MeV/;i) ^^^^ 
A„(w = 0) ° ' Ap(w = 0) 

For both tlie i/„ = 1, SD bands ^^'^Eg{l) and i93Hg(2b) 

An{uj = 0.50MeV/h) 



An{uj = 0) 

For the z/„ = 2, SD band ^^^Hg(2) 

A„(a; = O.SOMeV/^i) 
An(t^ = 0) 

IV. SUMMARY 



6% . (20) 



90% . (21) 



The u- and //-dependence of the nuclear pairing gaps of multiquasiparticle bands in 
well-deformed and SD nuclei are calculated under the PNC formalism, in which the block- 
ing effects on pairing are exactly taken into account. PNC calculations show that the 
cj-dependence of pairing gaps A for the u = (qp-vacuum) bands is weaker than what 
predicted in the particle-number projected HFB formalism. For the low-lying excited z/ > 2 
(~multiquasiparticle) bands, Ap's and A.„'s keep almost w-independent. As a function of 
seniority u, the bandhead pairing gaps A{u = 0, i/), decrease slowly with increasing u. 
Even for the highest seniority bands identified so far, the pairing gaps Ap{uj = 0, i/) and 
An{u! = 0, z/) remains larger than 70% of the bandhead value of the qp-vacuum band. 
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